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The validity of the pn-QRPA and -RQRPA descriptions of double 
beta decay transition amplitudes is analyzed by using an exactly solv- 
able model. It is shown that the collapse of the QRPA is physically 
meaningful and that it is associated with the appearance of a state 
with zero energy in the spectrum. It is shown that in the RQRPA this 
particular feature is not present and that this approach leads to finite 
but otherwise spurious results for the double beta decay transition 
amplitudes near the point of collapse. 
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The neutrinoless double beta decay (PPou) violates lepton number conser- 
vation and requires the existence of massive Majorana neutrinos Due to 
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this fact the detection of this decay has attracted considerable experimental 
effort @. The two neutrino mode of the double beta decay (Pfav), on the 
other hand, can be described as a second order process in the standard elec- 
troweak theory, its decay rate being independent of any new physics beyond 
the standard model. Both to predict and to analyze the data on double beta 
decay require a precise calculation of the various nuclear matrix elements 
needed to compute the corresponding half-lives. 

The quasiparticle random phase approximation (QRPA), including a par- 
ticle - particle channel in the residual interaction, can reproduce the experi- 
mentally determined two- neutrino double-beta decay {fifiiv) half-lives. This 
is so because the transition amplitudes are strongly suppressed for certain 
values of the force parameters. This result was tested by comparing the re- 
sults of the QRPA against an exactly solvable model || f| and including 
realistic residual interactions 

However, the predictive power of the QRPA is strongly diminished be- 
cause the ground state to ground state fifau transition amplitudes are ex- 
tremely dependent upon the structure of the adopted proton-neutron inter- 
action. For some critical values of the model parameters the QRPA collapses, 
i.e; the energy of the first excited QRPA state vanishes. It makes the theory 
nearly useless for some particular cases, the worst of which is the (Ifav decay 
of 100 Mo |, | . 

A renormalized version of the QRPA (RQRPA) JlO, llfl , which includes 
some corrections beyond the quasiboson approximation, has been reformu- 
lated recently (T^j and applied to the fifiiv decay fI3 |. Contrary to the QRPA 
in the RQRPA there is no collapse for any set of values of the force's parame- 
ters. It was presented as a reliable tool and was applied to the Pfov decay of 
100 Mo [ 13| . Similar results were found with the inclusion of proton-neutron 



pairing interactions [[14] . 

There are some concerns about the formulation of the RQRPA. The most 
important one is that it makes use of the exact commutation relations be- 
tween phonons, adding one-quasiparticle scattering terms to the ordinary 
quasi-boson approximation, but it does not include similar terms in the 
hamiltonian and in the transition operators. In the present letter it will 
be shown that the collapse of the QRPA correlates with the presence of a 
state with zero energy which is not found in the RQRPA spectrum. This 
shortcoming is not only a failure of the RQRPA but it is shared by others 
higher order approximations. Spurious excitations are exactly separated in 
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the RPA but higher order approximations can introduce spurious components 



which finally dominate the results [15 



The model hamiltonian | |TE| , |TT| , consists of a single particle term, a 
pairing term for protons and neutrons and a schematic charge- dependent 
residual interaction including particle-hole and particle-particle channels. It 
has been shown that this interaction, when treated in the framework of the 
QRPA, produces similar results as those obtained by using a G-matrix con- 
structed from the OBEP Bonn potential in reproducing single- and double- 
beta decay matrix elements Hl8j [T9| , |20f . In this letter we will consider the 
single shell limit (j p = j n = j) and monopole term (J = 0) of the residual 
interaction. As we shall show later on this model, which is not intended 
to accurately reproduce actual data, does indeed display the qualitative fea- 
tures of a realistic pn-QRPA calculation. Excitation energies, single- and 
double-beta decay transition amplitudes and ground state correlations, in 
this model, depend on the particle-particle strength parameter in the same 
way as they do in more elaborated calculations with many single particle 
levels and with more realistic interactions. The present case corresponds to 
beta decay transitions of the Fermi type. 

The schematic hamiltonian reads 

H = ^2e p ala p -G p SlS p + ^2e n al l a n -G n SiS n + 2xf3~ -f3 + -2kP~ -P + , (1) 

p n 

with 



sl = e 44/ 2 > s i = e 44/2, 

p n 

P~ = Eij < A T ~\j > 4«i p ~ = Eij < i\r'\j > a\at 



(2) 



4 = 4 % being the particle creation operator, a p = (—iy p ~ mp Oj p _ mp its 
time reversal and r~ the isospin lowering operator (r _ |n >= \p >). The 
parameters x an d k play the role of the renormalization factors g p h and g pp 
introduced in the literature 0, |], |^, |^| . 

The hamiltonian ([!]) can be expressed in terms of generators of an SO (5) 
algebra [23|]. It can be reduced to an SO (5) and isospin scalar |23| if 



its parameters are selected as 

e P = e n , x = 0, G p = G n = 4k. (3) 
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This SO (5) scalar limit was used to study the "pairing plus monopole 
model" many years ago [25, 2^]. If x 7^ t ne isospin symmetry is broken 



in the particle-hole channel and for 4k ^ G p (G n ) this symmetry is broken 
in the particle-particle channel. The Hilbert space is constructed using the 
eigenstates of the isoscalar hamiltonian, which are labeled by the number of 
particles N, the isospin T and its projection T z . If G ^ 4k the eigenstates 
have definite N and T z but not good isospin since the Hamiltonian mixes 
states with different isospin T. 

We have selected N n > N p and a large j to mimic the realistic situation in 
medium and heavy nuclei. In this letter we will use the following parameters: 

j = 19/2, N = 20, 1 < T z < 5, 

e p = 0.7MeV, e n = O.OMeV, G p = G n = 0.2MeV, \ = 0.025Me\/ . 

(4) 

k will be used as a running parameter. 
Fig. 1 

Lowest energy + states of different nuclei are shown in Fig. 1, for G = 4k, 
in an energy vs. Z plot. States are labeled by (T, T z ). It is clear that a simple 
model with both like particles and pn- pairing interactions can reproduce 
the qualitative form of the mass parabola. Fermi transitions (j3~ = £~ ) 
are allowed between members of the same isospin multiplet. In our single 
shell example, lowest + states in each odd-odd nuclei (N-l, Z+l, A) are 
the isobaric analog states corresponding to the states of the even-even nuclei 
with (N, Z, A) nucleons. Ground states of the initial (N, Z, A) and final (N-2, 
Z+2, A) even-even nuclei have different isospin and the Fermi-double-beta 
decay is forbidden in the isoscalar limit. 

Fig. 2 

The results shown by a full line, in Fig. 2, represent the excitation energy 
E exc of the lowest + state in the odd-odd intermediate nucleus (T z = 3) with 
respect to the parent even-even nucleus (T z = 4), as a function of 4k/ G. It 
is clear that when 4k/ G ~ 1.4 proton-neutron correlation dominates over 
proton-proton and neutron-neutron pairing correlations and the excitation 
energy goes to zero. This behavior is not a surprise and is similar to that 



found in the case of the pairing plus quadrupole Hamiltonian |27|. In this 
case, if the quadrupole-quadrupole interaction is strong enough, the system 
becomes permanently deformed. This was shown long time ago using the 



"pairing plus monopole" model [25, 26] which is a solvable two-level model 
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obeying the algebra of SO (5). In this letter we have applied the same ideas 
to the proton-neutron problem. 

It must be stressed that, in realistic cases, odd-odd nuclei in the lower 
energy sector of the mass parabola cannot have negative excitation energies, 
as compared to the even-even nuclei with the same mass, because it will 
be in contradiction with the main evidence for pairing effects in medium 
and heavy nuclei. It would also suppress completely the double beta de- 
cay because the single beta decay would be allowed between even-even and 
odd-odd nuclei. For the cases where the hamiltonian ([!]) predicts negative 
excitation energies, i.e: the above mentioned overbinding of odd-odd nuclei, 
quadrupole-quadrupole interactions and permanent deformations have to be 
considered. 

Now we turn to the QRPA, the RQRPA and their validity. After perform- 
ing the Bogolyubov transformations, separately, for protons and neutrons, we 
have obtained the qp-hamiltonian 



H= (e p - X P )N P + (e n - X n )N n + X^A + A 2 (AUt + AA) 

-X 3 (A^B + B*A) - A 4 (A t 5 t + BA) + X 5 B^B + A 6 ( J B t J B t + BB) 

(5) 

being e p = e n = Gfl/2 the quasiparticle energies and X p , X n the chemical 



potentials. Introducing the quasiparticle creation operators a P , at [[HJ and 
Q = (2j + l)/2, we have defined: 



At 

A 3 
A 5 



j=o 

M=0 



St 



a p <g> a f , 



M=0 



i=p,n 



4Q 



x{u 2 p vl + v 2 p u 2 n ) - K{u 2 p u 2 n + v 2 p v 2 n ) , A 2 = 4fi(x + K)u p v p u n v. 



AVt(x + n)u n v n {u 2 p - v 2 ) , A 4 = 4fi(x + K)u p v p {u 2 n - v 2 ) 



4fi 



X(%K + v 2 p v 2 n ) - K{u 2 p v 2 n + v 2 u 2 n ) 



A fi 



-A 2 



(6) 

The operators A, A^ together with their counterparts for identical parti- 
cles and B, B^ , N p , N n are the generators of the SO (5) algebra [^T] . 

The linearized version of hamiltonian (J5|) is obtained by keeping only the 
first line of Eq.(5), i.e. by neglecting the scattering terms proportional to B 
and I?t. Its solutions were discussed in a previous paper |28| . 

Finding the eigenvalues and eigenvectors of hamiltonian (|5]) requires the 
same algebraic techniques involved in solving the original hamiltonian. But 
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the complexity of the problem increases severely, due to the fact that neither 
the quasiparticle number nor the quasiparticle isospin projection (or equiva- 
lently the number of proton and neutron quasiparticles) are good quantum 
numbers. It implies that the dimension of the basis increases in two orders of 
magnitude. Technical details and a number of different examples are given 



elsewhere [24 



The dashed line in Fig. 2 shows the dependence of the excitation energy 
E qp for hamiltonian (0). It reproduces the exact results reasonably up to 
4k/G ~ 1.4. From there on it goes to zero instead of taking negative values. 

The QRPA matrix is a 2 x 2 one, with sub-matrices Aquva = 2e + Ai 
and Bqtzva = 2A 2 . The eigen-energy is E QRPA = [(2e + A x ) 2 - 4A 2 ,] 1 / 2 . It 
is shown as a large-dot line in Fig. 2. It becomes an imaginary number 
if 2A2 > 2e + Ai- It means that for this limit the zero-boson component 
of ground state ceased to be dominant P8|. The collapse occurs very near 
the point where the exact excitation energies become negative. It is a very 
important point. The overestimation of the proton- neutron correlations in 
the QRPA mimics the more complicated physics found in the exact case. 
However it gives a clear signal about drastic changes in the correlations, in 
this region of the parameters governing the residual pn-interaction. 

In the renormalized QRPA the structure of the ground state is included 
explicitly fllT] . The RQRPA energy Erqrpa is always real. Its value must 



be obtained by solving simultaneously a set of non-linear equations |13], 
It is shown as a small-dot curve in Fig. 2. This figure strongly resemble Fig. 
1 of ref. []13] and Fig. 2 of ref. ||, where the energy of the first excited state 



is plotted against the particle-particle strength parameter g pp . The curves 
for the QRPA and the RQRPA are quite similar to those shown here. The 
new feature discussed in this letter is that the exact excitation energies are 
closer to the QRPA energies rather than to the renormalized ones. 
The amplitudes M 2u are evaluated as: 

M ^ < 0/I/TIOa >< OAiriOi > (7 s 
M " - V E X -E 1 + A (7) 

We have selected A = 0.5 MeV. The values of M 2u , corresponding to 
the exact solution, are shown in Fig. 3 (full line) as a function of 4k/ G. For 
the other approximations we have diagonalized the hamiltonian (|5|) for the 
initial (N p = 6, N n = 14, ground state |0; >) and final (N p = 8, N n = 12, 
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ground state |0/ >) ground-states of the participant even-even nuclei. The 
overlap between the two descriptions of the intermediate states belonging to 
the double-odd nucleus [J was included in the above formula (Eq.(7)). The 
sum runs over the allowed intermediate states corresponding to this model 
space 
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Fig. 3 

In all the cases the curve corresponding to M 2 „ is very similar to that 
found in realistic calculations |3|, |5|, || |13| , including its cancellation near 
the collapse of the QRPA description. The RQRPA extends this curve far 
beyond the value of k at which the QRPA collapses. But the validity of this 
result can be questioned because the RQRPA missed the vanishing of the 
excitation energy, as we have discussed before. 

To conclude, we have presented an exactly solvable model which resembles 
the main features of the realistic models used to describe the structure of the 
nuclei involved in double beta decay processes. We have compared the exact 
values for the excitation energies and double beta decay transition amplitudes 
with those obtained with the approximate qp-hamiltonian, the QRPA and 
renormalized RQRPA ones. We have shown that the collapse of the QRPA 
correlates with the presence of a state with zero excitation energy. This state 
is present in the exact solution of the model. It is found that the RQRPA 
solutions do not show the presence of such state. As a direct consequence of 
this fact this approximation gives finite but spurious results for the transition 
matrix elements, which unfortunately are not supported by the exact results. 
The presence of spurious states, introduced by the renormalization procedure, 
could be responsible for such an odd behavior. The effect of these spurious 
contributions upon the transition amplitudes and the mixing of orders, in 
the sense of the order classfication in powers of 1/fi, found in the RQRPA 
wave functions are discussed in detail somewhere else p4]| . 
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Figure Captions 

Fig. 1: Lowest energy states of the different nuclei, labeled by (T,T Z ), 
for G = 4k. Allowed Fermi transitions are denoted by arrows. 

Fig. 2: Excitation energies as a function of 4k/ G. 

Fig. 3: (3(3 2 u transition amplitudes M 2v vs. 4k/G. 
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